Edge state effects in junctions with graphene electrodes 
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We consider plane junctions with graphene electrodes, which are formed by a single-level system 
("molecule") placed between the edges of two single- layer graphene half planes. We calculate the 
edge Green functions of the electrodes and the corresponding lead self-energies for the molecular 
levels in the cases of semi-infinite single-layer electrodes with armchair and zigzag edges. We show 
two main effects: first, a peculiar energy-dependent level broadening, reflecting at low energies 
the linear energy dependence of the bulk density of states in graphene, and, second, the shift and 
splitting of the molecular level energy, especially pronounced in the case of the zigzag edges due to 
the influence of the edge states. These effects give rise to peculiar conductance features at finite 
bias and gate voltages. 

PACS numbers: 73.63.-b, 85.65.+h, 72.80.Vp 
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I. INTRODUCTION 

Single-molecule nanosystems are in the focus of ex- 
perimental and theoretical investigations in recent years. 
One branch of the basic research in this field is concen- 
trating on the question of the current through molecu- 
lar junctions with metallic or semiconductor electrodes. 
Many phenomena, typical for nanoscale transport, such 
as Coulomb blockade, Kondo effect, vibronic and pola- 
ronic effects, to name a few, have been observed and 
explained. A comprehensive bibliography can be found 
in Refs.— —. However, despite the experimental progress 
and the theoretical efforts, the understanding of the prop- 
erties of single molecules coupled to metal electrodes, es- 
pecially their transport properties, is far from being sat- 
isfactory. One of the issues is still the often poor repro- 
ducibility of the experimental results. Unfortunately, the 
structure and quality of molecular junctions are not com- 
pletely controlled yet. One of the main problems is the 
size mismatch between metal electrodes and molecules 
and the impossibility to control the metal-to-molecule 
interface at nanoscale. Besides, gold, the most popular 
electrode material, has high atomic mobility and at room 
temperature the geometry of nanoelectrodes is not com- 
pletely stable. Thus, traditional electrodes show their 
principle limits and other materials for molecular elec- 
tronics should be considered. 

Carbon based materials, e.g., fullerene, nanotube or 
graphene, have the advantage of well controlled crystal 
structures, stability up to high temperatures due to sp2 
covalent bonds, and appropriate sizes. Production of ul- 
trathin epitaxial graphite films 4 opens new ways for elec- 
trode technology. 

Although the control of the atomic structure of 
graphene edges is also a challenge for present-day lithog- 
raphy, important progress has recently been achieved^ —. 



It is also important that carbon electrodes allow for many 
ways to anchor organic and inorganic molecules, thus be- 
ing promising for functional devices. Several types of 
carbon nanoelectrodes were suggested with different ge- 
ometry and dimensionality. One can mention, as an ex- 
ample, the study of electron transport across molecu- 
lar junctions with carbon nanotube electrodes (see, e.g., 
Refs. [IqIIiH and references therein). 

The advantage of plane graphene (single- or multi- 
layer) is that graphene sheets are suitable for nanoscale 
lithography 4 ^— and can be considered as a base for molec- 
ular integrated circuits. Thus we focus on molecu- 
lar junctions with plane graphene electrodes (Fig.[T]), 
where a molecule is coupled to the graphene edges (here 
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Figure 1: (Color online) Example of a plane graphene molec- 
ular junction (upper panel) and its schematic representation 
(lower panel) of the considered single-level model for zigzag 
edge. 
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a variety of molecules and anchoring groups can be 
used). The theoretical investigation of such structures 
was started recently for all-carbon junctions^ 2 -"— and 
junctions with different organic molecules.— In partic- 
ular, one can mention the experimental^ and theoreti- 
cal & 14 i 16 — investigations of linear atomic carbon chains 
between single-layer graphene electrodes and first exper- 
imental observation of the current through the molecular 
junction with few-layer electrodes. 25 

The other advantage of graphene electrodes is that the 
molecular gating problem could be solved. Indeed, large 
metal electrodes screen the external gate potential and 
makes it almost impossible to shift molecular levels in a 
controlled way. Oppositely, in plane structures with thin 
graphene electrodes the gate potential can be used rather 
efficiently. 

The simplest possible case of a molecular bridge with 
graphene electrodes is a junction formed in a nanogap 
with armchair or zigzag edges. The zigzag edge case is 
shown schematically in Fig.[TJ The armchair edge is ob- 
tained if one cuts a graphene sheet in the perpendicular 
direction. The transport properties of such junctions are 
determined by the peculiarities of the edge Green func- 
tions of semi-infinite graphene electrodes. In the case 
of zigzag edge the main features are influenced by edge 
states. In this paper we address the question of how these 
properties of graphene electrodes affect charge transport 
through the molecular bridge. Hence we do not focus on 
molecular specific properties and replace the molecule by 
one spin-degenerate noninteracting level. The electrodes 
are described in the framework of the 7r-electron tight- 
binding model. 26 This approach is complementary to an 
ab initio one and allows to obtain physically transparent 
results. 

Below we consider the tight-binding model of graphene 
molecular junctions with armchair and zigzag edges 
(Sec. HI). We calculate the edge Green functions and lead 
self-energies in Sec.lIIII Then, in Sec.|lVlthe edge depen- 
dent level shift, splitting and broadening are considered. 
Finally, in Sec. [V] we discuss the current at finite voltage 
and the gating effect, and give conclusions in Sec.lVTl 

II. THE TIGHT-BINDING MODEL 

We write the tight-binding Hamiltonian of our system 
(for fixed spin a =t or a the spin index is omitted) 

as 

H =(e + eVg + e<p )rfd + {V£c\ n d + V^c sn ) 

s=L,R 

+ ^2( eV 9l + e( Ps) C ] si C si + ^2Uj C ] si C sj . (1) 

s; i s; ij 

Here d) , d are the creation and annihilation operators 
for a molecular level, while the c^, c S i operators describe 
the lattice tight-binding local orbital at the i-th site in 
the s-th (left or right) electrode. The molecular level 



is assumed to be coupled only to one edge site of the 
graphene lead with site index n (Fig. [I]) and Vb is the cou- 
pling matrix element. We denote the tight-binding hop- 
ping matrix elements tij between lattice sites by Uj = t 
for nearest neighbor sites, and = t' for next-nearest 
neighbor sites of the lattice. V g is the gate voltage act- 
ing on the molecular level and V g \ is the gate voltage, 
applied to the leads. At finite bias voltage V (defined 
by the left and right electrical potentials, V = pl — Pr) 
the energy of the molecular level is shifted. In the linear 
approximation this shift is described by a parameter 77: 
^0 = <fR + vi^L — Vr)-, where < 77 < 1 characterizes 
the symmetry of the voltage drop across the junction, 
and r] = 0.5 stands for the symmetric case (when it is 
convenient to use pl = V/2 and pr = —V/2). The gen- 
eralization of this model to the case of many molecular 
levels, as well as interacting molecular levels, is straight- 
forward. 

Note that while the lead gate potential V g i and the 
bias potentials p s enter the Hamiltonian (pQ) on an equal 
footing, the physical sense and the effect of these two po- 
tentials are different. The gate potential shifts the energy 
levels in the finite-size region of the electrodes near the 
molecule, where it is applied, but does not change the 
Fermi level, which is determined by the large equilibrium 
electrodes outside this region. The bias potential shifts 
additionally the energy distribution of electrons in the 
electrodes according to 



s W exp((e-e<p a )/k B T) + l v 7 

with temperature T. 

The gate potentials V g i and V g of the leads and the 
molecule can generally differ, but in this paper we con- 
sider the case V g = V g i being a good approximation for 
plane structures. 

III. EDGE GREEN FUNCTIONS AND 
SELF-ENERGIES 

There are two electronic properties that determine the 
main peculiarities of graphene edges as electrodes. First, 
the energy dependence of the local density of states in 
bulk graphene p(e) has a minimum near the Fermi en- 
ergy of undoped graphene (to be chosen as the zero en- 
ergy) and is almost linear near this point p(e) oc |e|. The 
coupling of molecular electronic levels to the electrodes 
is characterized by the lead self-energy E(e), which is 
proportional to the edge Green function of the electrode 
G nm (e). If a single molecular level is coupled to one 
edge atom of the electrodes, as in our case, the retarded 
(lesser) self-energies are 

Sf«)(e) = |^o| 2 G^<)(e), (3) 

where G s ,nn is the Green function of the 5-th lead at the 
site n, to which the molecular level is connected. 
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Figure 2: (Color online) Edge Green function G nn (e) for the 
armchair edge. 



The self-energy is a complex function, with 
the imaginary part known as the level-width 
function T(e) ex — ImE jR (e) ex |Fo| 2 p( e ) (actually 
Pn(c) = —^ImGnnic) at lattice site n). The level- 
width function determines, in particular, the broadening 
of molecular levels and the maximal (resonant) current 
through these levels. Thus, the current will strongly 
depend on the relative position of the molecular level, 
the Fermi level and the Dirac point (the energy of 
minimal p(e)). This is true for both armchair and zigzag 
edges. Moreover the zigzag edge supports so-called 
edge states j£r2!L which are localized near the edge with 
energies near the Dirac point. 

We performed numerical calculations of the edge Green 
functions and lead self-energies (j3j) for armchair and 
zigzag graphene edges using the iterative metho d 31 ! 32 and 
the eigendecomposition based method^ 3 - The results for 
the edge retarded Green function are presented in Figs. 
[2] and G3 We show here only the diagonal part G^ n of 
the full matrix Green function, required to calculate the 
self-energy. We checked the influence of additional next- 
nearest neighbor coupling t' (Fig. [3b) and found that the 
main presented results are qualitatively the same. Thus 
we focus on the case t' = in this paper. 

In both, of the armchair and zigzag cases we found the 
linear energy dependence of the imaginary part (density 
of states) at small energies, which reflects the bulk prop- 
erties of graphene. In the zigzag case an additional fea- 
ture in the density of states is present near zero energy, 
and the real part of the self-energy is quite different from 
the armchair case. Because of the weak energy depen- 
dence on the quasi-momentum (along the edge) the edge 
states form a flat energy band with a ^-function type den- 
sity of states giving rise to the divergence close to e = 0. 
Correspondingly, the real parts of G nn (e) at the edge and 
E(e) have a singularity at this energy (Fig. [3]). Since the 
real part of the self-energy renormalizes the energy of the 
molecular levels, e' ~ eo + ReS(eo), the energies of levels 



at low energies are strongly shifted. Moreover, the spec- 
tral function, and hence the dressed energy levels, can 
be split because the sign of ReE(e) changes close to the 
singularity point. 



IV. LEVEL SHIFT, SPLITTING AND 
BROADENING 



Using the calculated self-energies we investigated the 
spectral function 



A(e) = -2lmG R (e) 



(4) 



of a single level, coupled to armchair and zigzag elec- 
trodes. The retarded Green function of the level is 



G R (e) = 



6-6 -£2(6)-£«(6)' 



(5) 



where the lead self-energies S^(e) and S^(e) are calcu- 
lated from the Green functions of the leads using expres- 
sion (j3]). 
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Figure 3: (Color online) Edge Green function C7 n n(e) for the 
zigzag edge with (a) only nearest neighbor hopping t and (b) 
additional next-nearest neighbor hopping t' — O.lt. 
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The renormalized resonant level position e' is de- 
termined from the maximum of the spectral function, 
Eqs.gEJ): 



e' - e - Re [Sf (e') 



0. 



(6) 



For the armchair case (Fig. [2]) the graphical solution of 
this equation has one solution in most cases, but three 
solutions at eo ~ d=t. In the zigzag case (Fig. [3]) there are 
always two solutions because of the sign-changing singu- 
larity at small energies. The existence of two solutions 
means that the original level is split into two sublevels. 
This splitting can be understood as a result of the hy- 
bridization of the molecular level with the edge states. 
The magnitude of the splitting (the distance between the 
levels), as well as the shift and broadening, are controlled 
by the coupling Vb of the molecular level to the leads. Be- 
low we consider mainly the most interesting case Vb <C t 
and choose Vb ~ O.lt, which is typical for covalent cou- 
pling of organic molecules to the leads. 

The broadening of the resonant level with energy e' is 
determined by the imaginary parts of the self-energies, 



r = -2Im [E?(e') + £g(e' 



(7) 



It is mainly given by the graphene bulk density of states 
(with some quantitative variations in the armchair and 
zigzag cases). Additionally there is a sharp peak of the 
density of states at the energy of the edge state; however, 
the molecular levels are shifted from these energies. 

Figures |4] and [5] show A(e) for armchair and zigzag 
lead termination, respectively, as a function of the en- 
ergy for different original level positions eo- For com- 
parison we present also the curves in the wide-band limit 
(E^(e) = — (i/2)T s = const), commonly chosen for metal 
leads. As we already anticipated, the main results are a 
energy level position dependent broadening for both edge 
types and a level shift and splitting mainly for the zigzag 
edge. Of course, the results also depend on the coupling 
V . 

In the armchair case (Fig. [4]) the main effect on A(e) at 
small level energies (eo < t) is the energy-dependent level 
broadening due to the imaginary part of the self-energy, 
but the level shift due to the real part of the self-energy 
is also clearly observed (Fig.HJ)). The peak at eo = 
is very sharp, because the imaginary part vanishes at 
e = 0. At larger level energies the spectral function has 
a more complex form (Fig. 0k) • The levels at eo ~ ±£ 
are split, because ReS R changes sign at e = ±t (Fig. [2]). 
However due to the maximal level broadening at these 
energies, the splitting is not clearly pronounced and the 
levels overlap. At larger energies the broadening and shift 
decreases. 

In the zigzag case (Fig. [5} the singularity of ReS R at 
zero energy changes the picture drastically, especially for 
small level energies. First of all, the spectral density is 
"repelled" from small energies. If the level is originally 
at eo = 0, it is split into two sublevels. At large |eo| the 
small part of the integral spectral density of the original 



level is split and the second sublevel appears at the other 
side of the e = point. For large energies of the original 
level, the second sublevel is close to zero energy (position 
of the edge state), but its spectral weight is small. 

We summarize these findings in the 2D plot, Fig. [6] 
Here the spectral function of a single level is shown color- 
coded as a function of the energy e and the level position 
eo in three cases: the wide-band limit (for comparison), 
the armchair leads and the zigzag leads. All the three 
effects: level broadening, shift and splitting are visible in 
Fig.[6fb) and Fig.[6fc). The strong energy dependence of 
the edge Green functions and correspondingly of the lead 
self-energies of graphene electrodes leads to the qualita- 
tively different behavior of the molecular levels, coupled 
to such electrodes, compared to the case of wide-band 
electrodes. The shift of the level position and the broad- 
ening depend strongly on the energy of the unperturbed 
molecular level eo- Besides, in the zigzag case, the edge 
states of the leads are hybridized with the molecular level, 
and the molecular level is split into two sublevels. 

Such unusual interplay of level shift and broadening 
is expected to result in peculiar transport properties 
through graphene molecular junctions, compared to the 
case of metal wide-band electrodes. 
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Figure 4: (Color online) Spectral function of a single level 
coupled to armchair edges. The equidistant dashed lines show 
the position of the level in the wide-band limit with Vo = O.lt. 
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V. CURRENT AND DIFFERENTIAL 
CONDUCTANCE: GATING AND EDGE STATE 
EFFECTS 



Finally we calculate the current and differential con- 
ductance with the assumption that graphene electrodes 
(including the edge states) are kept in equilibrium, but 
at different electrical potentials due to the bias voltage 
V. This condition can be easily fulfilled if coupling of a 
molecule to the leads is weak enough. We follow the for- 
mulation pioneered by Meir, Wingreen and Jauho*24~— 
The current from the left (s = L) or right (s = R) lead 
into the central system is described by the expression 
(here we consider the case of spin-unpolarized leads) 

j -=!/! Tr f r 'M<c < H+ 

[G*(e)-G A (c)])}. (8) 

Here fg(e) is the equilibrium Fermi distribution func- 
tion in the 5-th lead, Eq. ((2]), T s (e) = -2Im£f (e) is the 
level- width function, G R (e) is the retarded function of 
the level, as defined by Eq. ©, G A (e) = [G R (e)]*, and 
G < (e) is the lesser Green function. It can be found from 
the Dyson-Keldysh equation in the integral form, 



G<(e) = G R (e) (£<(e) + £<(e)) G A (e), 



(9) 
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Figure 5: (Color online) Spectral function of a single level 
coupled to zigzag edges. The dashed lines show the position 
of the levels in the wide-band limit with Vb = 0.1£. 
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Figure 6: (Color online) Color-coded strength of spectral 
function of a single level coupled to wide-band leads (a), arm- 
chair leads (b) and zigzag leads (c) . The scale of the spectral 
function is shown at the right. 



where the lesser self-energy of the noninteracting leads is 
E<(e) = ;r s (e)/ S °(e). (10) 

The current is a function of the bias voltage 
V = (fL — PR- The second parameter is the gate volt- 
age V g i which is assumed to shift both the energy level 
e = eo + eVgi of the molecule and the energy levels in 
the graphene electrodes, which results in the shift of the 
self-energy E(e) = E(e - eV g i). 

In Fig. [7] the differential conductance dl/dV for eo = 
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Figure 7: (Color online) Differential conductance dl/dV as a 
function of bias and gate voltage for a single level at eo = 
coupled to wide-band leads (a), armchair leads (b) and zigzag 
leads (c), with V = O.lt and k B T = O.Olt. The scale of the 
differential conductance is shown at the right. 



is shown as a function of bias and gate voltage for arm- 
chair and zigzag leads in comparison with the wide-band 
limit. The main features of conductance are explained 
by the level modification properties, considered previ- 
ously. For example, at small V and V g i (in the center) 
the differential conductance has a minimum rather than 
a maximum, contrary to the metal wide-band case with 




d) 



Figure 8: The energy level diagrams corresponding to the 
points a)-d) marked in Fig.JTJc) of the differential conductance 
for zigzag electrodes. The vertical axis corresponds to the 
energies. The linear energy dependence of the bulk density of 
states near the Dirac point is shown schematically for the left 
and right leads (the singular edge-state density of states at 
the Dirac point is not shown). See the text for explanations. 



energy-independent lead self-energy. 

In the case of armchair electrodes the level is not 
shifted and the small conductance at small bias and gate 
voltage is mainly explained by the small T(e = 0) that 
gives also the small current I ocTV. If we apply the gate 
voltage, the maximum of differential conductance occurs 
at V = ±2V g i, when the Fermi energy crosses the level 
position, and the current at this point is determined by 
T(e = Vgi), which is larger at larger V g \. On the other 
hand, the level broadening is also larger in this case. The 
regions of negative differential conductance at large gate 
voltage and large bias voltage appear because of the de- 
creasing density of states at e > t. 

In the case of zigzag electrodes, this picture is compli- 
cated by the edge-state effects, removing levels at small 
energies and suppressing the transport at small voltages. 
The additional structures appear due to the existence of 
the split sublevel. It can lead to interesting effects, when 
the current at small voltages is determined by the molec- 
ular levels with large energies, while the small energy 
levels are shifted to larger energies. Besides, the level 
broadening and shift take place for larger |eo| similar to 
the case of the armchair edge. For a further analysis of 
the dl/dV characteristic in the zigzag case, let us con- 
sider the energy level diagrams (Fig.|8j) corresponding to 
the points a), b), c), d) on Fig.[7h. At point a) the level 
is split into two equivalent sublevels at finite energy, so 
that the linear conductance is small. At finite gate volt- 
age (point b)) the Fermi level in the leads is in resonance 
with one of the sublevels and a distinct local conductance 
maximum is observed at zero voltage. At zero gate volt- 
age, but finite bias voltage (point c)), the Dirac points 
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and the energies of the edge states are shifted away from 
zero energy. As a result, the original energy level with 
eo = is not shifted, but produces two additional side 
levels at positive and negative energies. The current flows 
mainly through the unshifted level and is a monotonous 
function of the voltage, because it is proportional to the 
density of filled states from one side and the density of 
empty states from the other side. Similar behavior is ob- 
served at this voltages in the armchair case, because edge 
effects do not play an essential role. Finally, at point d), 
with \V\ ~ |2V^|, the Fermi level in one of the leads 
crosses the energy level and produces the resonance line. 
This is similar to the wide-band case, but the current and 
conductance are larger at larger voltages, because of the 
energy dependent coupling to the leads. 

VI. CONCLUSIONS 

In conclusion, we showed that gate effects in graphene 
molecular junctions are determined by two main effects: 
the strong energy dependence of the lead self-energies 
and the shift of these energy-dependent self-energies rel- 
ative to the Fermi level, together with the usual shift of 
the molecular levels. In particular, in the zigzag case we 



demonstrate the splitting of the molecular resonant state. 
In the biased case with bipolar graphene leads discussed 
above, removing the nanogap and the molecular level re- 
covers a graphene pn junction, where the transmission 
becomes reflectionless at normal incidence. This perfect 
transmission, decaying with the increase of the incidence 
angle, is known as Klein tunneling, which has been pio- 
neered for graphene by spin-independent studie o 37 i 38 and 
recently generalized to spin-dependent cases j 39 i 4Q How 
the picture of Klein tunneling would be modified by the 
presence of the nanogap with a (single) bridging molecu- 
lar level is naturally an interesting question that deserves 
future investigation. 
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